Blurred femtoscopy in two-proton decay 
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We study the effects of final state interactions in two-proton emission by nuclei. Our approach 
is based on the solution the time-dependent Schrodinger equation. We show that the final relative 
energy between the protons is substantially influenced by the final state interactions. We also show 
that alternative correlation functions can be constructed showing large sensitivity to the spin of the 
diproton system. 

PACS numbers: 23.50.+Z, 27.20.+n, 27.40.+Z, 29.40.Cs 



Two-proton emission has been observed for numerous 
excited states in nucleLpopulated both in (3 decay and in 
nuclear reactions 

[E H S H II • Although these decays are 
thought to be sequential one-proton emissions proceeding 
through states in the intermediate nucleus [1, 0] , there is 
an intriguing possibility [|[ that the diproton ( 2 He) cor- 
relation may play an important role in the mechanism of 
the two-proton emission. This has been nicely demon- 
strated in the analysis of the two-proton decay of the 
21+ isomeric state in 94 Ag, where 19 events were clearly 
assigned to the simultaneous emission of two correlated 
protons Q. The traditional idea of diproton radioac- 
tivity is due to the pairing effect. Two protons form a 
quasiparticle (diproton) under the Coulomb barrier and 
this facilitates penetration. In a more formal description, 
one has a system with two valence protons in the same 
shell and coupled to J v = + . This question, being still 
open, continues to motivate studies in this field. In order 
to assess this information, it is necessary to understand 
final state interactions between the protons and between 
each proton and the daughter nucleus. 

In nuclear decays the emission of correlated, identical, 
particles is sensitive to the geometry of the system. Mea- 
surements of correlation functions are often performed 
with charged particle pairs, which interact via the short- 
range nuclear interaction and the long-range Coulomb 
interaction and they also interact with the remaining 
source. As a result, theoretical corrections are needed to 
subtract the final state interactions (FSI) before one can 
extract any useful information about the emitting source 
from the measurements 

@, [13, EH E2 - At first sight, the 
FSI can be regarded as a contamination of "pure" par- 
ticle correlations. However, it should be noted that the 
FSI depend on the structure of the emitting source and 
thus provide information about source dynamics as well. 

Two-proton decay in s-wave states can also be used for 
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testing quantum mechanics versus local realism by means 
of Bell's inequalities [HI]. Since the final state of the two 
protons can be either in a singlet or in a triplet state, 
their wavefunction is spin entangled. The identification 
of the spins of the proton in two detectors separated 
far away would be useful to test the Einstein-Podolski- 
Rosen (EPR) paradox [30]. In fact, these tests should 
be performed in different and complementary branches of 
physics to avoid the loopholes encountered in photon ex- 
periments. The advantage of using massive Fermions to 
test Bell-type inequalities is that the particles are well 
localized and the spin state of the pair can be well estab- 
lished by measuring the internal energy of the two-proton 
system. However, the validity of this method highly de- 
pends on our ability to treat FSI. Coincidence measure- 
ments of the two proton momenta require knowledge of 
FSI in order to extract information about their original 
wavefunction. Here will propose a new method to cal- 
culate FSI based on the numerical solution of the time- 
dependent Schrodinger equation. We hope with that to 
get a quantitative estimate of the FSI and how they can 
be used to address the points raised above. Though 
tests of Bell's inequality using proton-proton spin cor- 
relation in low energy scattering conformed to quantum 
mechanics (see, e.g., [Hj]), one would still like to find 
other means to perform the verification of the complete 
nature of quantum mechanics. 

One should distinguish this work from Hanbury- 
Brown-Twiss (HBT) studies in high enery nucleus nu- 
cleus collisions. Indeed, in the case of HBT, the whole 
game is played by FSI. FSI are the means by which one 
can determine information about the source. In this case 
FSI is not viewed as a " contamination" . The impor- 
tant difference in our case is that we are interested in a 
case where the protons are not emitted chaotically like in 
the case of HIC (Heavy-Ion Collisions). We are studying 
nuclear structure. In HIC (and HBT) protons are as- 
sumed to be emitted independently and chaotically from 
the source (any information about their initial spins is 
lost and they are assumed to be "evaporated" with a boil- 
ing pot). Then, there are no initial state correlations and 
FSI make the whole physics. In the case of two-proton ra- 
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dioactivity, the emission of the two protons is not chaotic 
because their correlation function keeps memory of their 
spin admixture and wave function in the parent nucleus. 

We consider first a single proton described at the initial 
time by a localized wave-packet -0o( r i)- The probability 
amplitude to find the proton at the detector with mo- 
mentum pi is given by 

A(p 1 ,r 1 ) = [dr X i+) {Pi,r)K(r,ri)Mn), (1) 

where X^Hpij 1 *) i s an asymptotic outgoing Coulomb 
wave with energy E = p 2 /2m p , and K (r, r-J is the prop- 
agator which accounts for the time evolution of the par- 
ticle from the source to the detector. 

We now look at the case of two-protons interacting 
with the residual nucleus and between themselves. We 
will consider the distortion caused by the Coulomb plus 
nuclear interaction between each proton i with the nu- 
cleus, Vc(ri) + VAf(ri), and between themselves, v}?(r) + 
i^(r), where r-j is the coordinate of proton i, and r is 
their relative coordinate. The proton-nucleus interaction, 
VAr(ri) yields smaller final state interaction effects than 
the Coulomb counterpart. 

We adopt a classical description of the center-of-mass 
motion for the two-protons and solve the time-dependent 
Schrodinger equation for the relative motion between 
them. The Coulomb field that distorts the relative mo- 
tion of the particles is given by 



V c {t) = Ze 2 



1 



1 



|n -R(t)| |r 2 -R(i)| R(t) 



, (2) 



where Z is the charge of the daughter nucleus and ri 
and r 2 are the positions of the protons with respect to 
the center of the nucleus (nuclear recoil is neglected). 
Vc(t) acts on the relative position r = r 2 — i"i through 
the transformations ri = R r/2 and r 2 = R + r/2. 

One can perform a multipole expansion of this inter- 
action and for r smaller than R(t) one can express the 
result in terms of a multipolc-dependent effective charge, 

= e (— 1/2) L + (1/2) L where L is the multipole de- 
gree. The dipole field (L = 1) is only important for 
particles with different charge-to-mass ratios while the 
quadrupole field is dominant when these ratios are equal 
(e.g. for two- proton emission) . For the quadrupole inter- 
action, eL=2 = e/2 and 



Vc(t) 



Ze 2 



2 R 3 (t) 



P 2 (cos ( 



(3) 



where 9 is the angle between R and r and P 2 s the Leg- 
endre polynomial of order 2. We now assume that the 
protons are produced simultaneously and nearly at rest 
at position 2ao and time t = 0. Their center-of-mass 
follows a radial trajectory described by 

R(t) = y (coshw + 1) , t = ^ (smhw + w) , (4) 



where the asymptotic velocity is given by v — y/E/m p , 
E is the two-proton decay energy, and a = e 2 /2E. This 
assumes that the relative energy between the protons is 
much smaller than E, which is not a good approximation, 
as we will show later. It is important to notice that eqs. 
E] only account for the motion of the protons after they 
emerge from inside the nucleus through the Coulomb bar- 
rier and propagate from the closest distance 2ao to infin- 
ity. Hence, our calculations neglect what happens during 
the tunneling process and treat only the external mo- 
tion. Hence, neglecting the proton-nucleus strong FSI is 
justified. 

We can still use eq. [T] to calculate the probabilities 
for relative motion of the protons, with the wavefunction 
for the relative motion given by \&(r) = K (r, r ) i^q(yq). 
In the time dependent description, at time t this wave 
function can be expanded in spherical harmonics 



* (r) = - ui m (r, t) Yi m (r) , 



(5) 



and the Schrodinger equation, describing the time evolu- 
tion of the relative motion between the protons can be 
solved by the finite difference method, calculating the 
wavefunction at time t + At in terms of the wavefunction 
at time t, according to the algorithm 
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where t = hAt/m p (ArY 
tor is defined as 



The second difference opera- 



A (2 M;2(t) = u^\t) + u^\t) 2ttj2(t), (7) 

with (t) = ui m {vj , t) , where rj is a position in the 

y {r.j) is the 



radial lattice. In eq. [51 U — v}? (rj) 
Coulomb+nuclear interaction between the two protons as 
a function of their distance, rj, and the function Si> m '-i m 
is given by 

Si> m >;im{r,t) = V* (Y Vm , \V C {r , t)\Y lm ) u Vm , (r,t) . (8) 



This method of solving the time-dependent equation is 
the same as used in ref. [17J for studying reacceleration 
effects in breakup reactions in nucleus- nucleus collisions 
at intermediate energies. A grid adequate for our pur- 
poses has 5000 spatial mesh points separated by 0.1 fm 
and 2000 time mesh points separated by 0.5 fm/c. 

We use the quantization-axis along the R(t) ccntcr-of- 
mass radial trajectory. As a consequence, P 2 (cos 9) = 
^/47r/5y 20 {9, (f>), and one only needs to consider the 
m = component of the spherical harmonics implicitly 
contained in the potential Vc(t). The initial I = state 
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cannot develop a final I = 1 component, and only I = 
(s-waves) and / = 2 (d-waves) will be present in the fi- 
nal state. Higher I values will be small and need not be 
considered. 

The proton-proton potential is taken as (r) + 
V C? ( r ) = e 2 / r + vo(b/r) exp (— r/b). The set of param- 
eters v Q = -46.124 MeV and b = 1.1809 fm yields the 



proton-proton scattering length, a v 



-7.8196 fm and 



the effective range po = 2.790 fm, in accordance with 
experimental data. But we choose a higher absolute 
value of Vq which allows the presence of a single weakly 
bound s-wave state. We use this localized wavefunction 
for the relative motion of the two protons in the initial 
state: uq = ui=q (r,t = 0). This is an artifact of the 
numerical method chosen as to allow for a localization 
of the initial wavefunction. The observables associated 
with the final state will depend on the binding energy, 
reflecting the dependence on the initial average separa- 
tion between them. The average initial separation, ro, 
and the binding energy, B, are approximately related by 
r = h(ABm p y 1/2 . 

As time evolves the initial state will acquire compo- 
nents in the continuum due to the action of the inter- 
action Vc(t). The continuum component propagates as 
a wavepacket which moves away from the source with 
a final asymptotic momentum p. The continuum wave- 
function is obtained by removing the bound-state part 
from the solution of eq. [5] 

* c (*)=7V[*(t)-<*(*)|* ) *o], 

where Af normalizes the continuum wavefunction, f c , to 
unity. 

The probability amplitude to find the protons with a 
final relative momentum p is given by 

*4(P) = (X (+ W)|*(r,<^» 



where 

Hi (pr) = exp 



ui=2(r,t - 
In 



■ oo)Ho(pr)dr 
oo)H 2 (pr)dr, (9) 



i( F +y-i ) ln(2pr) + a l ] (10) 



is the asymptotic Coulomb wavefunction for angular mo- 
mentum I, with p=hk, rj — e 2 /hv r , v r — y/ 'AE r /m p is 
the asymptotic relative velocity, and E r = h 2 k 2 /m p the 
relative energy. 

We consider two-proton decay from 45 Fe with a de- 
cay energy of 1.1 MeV. In figure [T] we show the re- 
sults for the average relative motion energy (E r ) = 
(A \p 2 /m p \ A) I (A\A) as a function of the average initial 
distance between the two protons. The dashed curve is 
the final proton relative energy if their mutual interaction 
is neglected. The dotted curve includes the Coulomb re- 
pulsion between the protons and the solid curve includes 
both their Coulomb and nuclear interaction. 
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FIG. 1: Average relative motion energy of the two protons 
as a function of the average initial distance between them. 
The dashed curve is the final proton-proton relative energy if 
their interaction is neglected. The dotted curve includes the 
Coulomb repulsion between them and the solid curve includes 
their Coulomb and nuclear interaction. 



The physical reasons for the results shown in figure [T] 
are transparent: two strongly interacting particles emit- 
ted from the volume of ro ~ 2 fm, when their "own mean 
radius", reflected by their attractive strong potential vn, 
is about b ~ 1 fm, should sufficiently "feel" each other 
through mutual attraction. On the other hand, when 
the emission zone is much larger than the particle mean 
radius, for instance if the radius is ro ~ 6 fm, the con- 
tribution of the short-range attraction is negligible and 
only the long-range Coulomb repulsion acts, as is seen 
in fig. [T] It is also clear that the final relative energy 
increases when the distance ro decreases. But as ro in- 
creases the contribution of the "tidal" Coulomb interac- 
tion of the diproton with the daughter nucleus decreases 
much faster than the contribution of their own mutual 
Coulomb repulsion. 

It is also important to notice that all the contributions 
to the FSI energies are not small compared to the decay 
energy E. This is contrary to our initial assumption used 
to justify our dynamical model. Hence, the results point 
to the important conclusion that final-state interactions 
are very important in determining the relation between 
the proton energies and the spatial distribution of the 
protons in the decay process. The FSI contributions due 
to the Coulomb tidal interaction depend on the square of 
the charge of the daughter nucleus while the contribution 
of the strong force between the protons is approximately 
independent of the nuclear mass. Thus, for lighter nu- 
clei (e.g. 18 Ne), the dashed curve in figure Q] becomes 
negligible. 

Further aspects of the spatial configuration of the pro- 
tons can be obtained in light of the usual discussion in 
terms of two-particle interferometry, or correlation func- 
tions. The relation between the center-of-mass coordi- 
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nates and the laboratory are given by 
(ri +r 2 ) 



R 



P = Pl +P2, 



r = ri - r 2 , 

_ (Pi - P2) 



The probability amplitude to find one of the protons with 
momentum pi is given by Ai (pi, ri) = A (pi — P/2, ri). 
The center of mass momentum, P, is set by the decaying 
energy and the assumption that it follows an outgoing 
radial motion, i.e. P = Ry E/m p , where R is the unit 
vector along the radial direction. 

Next we show that one can disentangle the contribu- 
tions of singlet and triplet spin final states of the two- 
proton system by measuring momentum correlations. 
This will prove to be a useful method since a direct mea- 
surement of the spin orientations of each proton is by 
far more complicated. The application of the method 
is very general as it only relies on measured quantities, 
independent of the models for the treatment of FSI. 

The protons are identical particles and their detec- 
tion requires the consideration of their quantum statis- 
tical properties. If proton 1 is detected with momen- 
tum pi and proton 2 is detected with momentum p 2 , 
the probability amplitude for this is given by product 
Ai (pi,ri)^42 (p2,r 2 ). Because of the indistinguishabil- 
ity of the particles, the probability amplitude must be 
symmetric with respect to the interchange of two par- 
ticles if they are in a spin-singlet state (S = 0), and 
antisymmetric if they are in a spin-triplet state (S = 1). 
The normalized probability amplitude becomes 

A (±) (pi,p 2 ,ri,r 2 ) = -j= [A! (pi,ri)^ 2 (P2,r 2 ) 

±Ai (p 2 ,ri)^2 (Pi,r 2 )] , 

where the plus or minus sign is the spin-singlet and spin- 
triplet state, respectively. 

The two-particle momentum distribution P(pi,p 2 ) is 
the probability to measure a nucleon having momentum 
Pi in coincidence with the measurement of the other nu- 
cleon having momentum P2. It is defined as 

P(pi,p 2 ) = Jd 3 nd 3 r 2 A<+' (pi,p 2 ,ri,r 2 ) 



±MA { >(pi,p 2 ,n,r 2 ) 2 



(11) 



where M. is the mixing parameter, determining the rel- 
ative contribution of the triplet state. The correlation 
function C(pi,p 2 ) is defined as the ratio of the prob- 
ability for the coincidence of pi and P2 relative to the 
probability of observing pi and P2 separately, 



C(pi,p 2 ) 



P( P i,p 2 ) 
Pi(pi)P 2 (p 2 )' 



Let us assume for the moment the that protons sud- 
denly emerge from the nucleus and that their intrinsic 



wavefunction is in a pure entangled state. If their wave- 
function is approximated by exp(ip 1 .r 1 ) exp(ip 2 .r 2 ) ± 
exp(zp 2 .ri) exp(ipi.r 2 ) it will lead to destructive or con- 
structive interferences. If we also assume a gaussian 



source of size rg = 



the correlation function with- 



out final state interactions would be given by 

C(pi, P2 ) = C(q)=f (Ap) [1 ± exp (-Ap 2 r 2 Q /h 2 )] , 

(12) 

where Ap = p = \pi — p 2 |/2. All other features of the 
reaction mechanism are included in the function / {Ap). 
One can approximately account for the Coulomb inter- 
action between the protons by using a Gamow function 
for/(Ap): 



f(Ap) 



exp (27rry) — 1 



(13) 



but no such simple estimate exists for the effect of the 
nuclear interaction. 

According to eq. for Apr^/h <C 1 one should 

be able to see a destructive interference for triplet fi- 
nal states and constructive interference for singlet final 
states. It is thus appropriate to redefine the correlation 
function in terms of the relative momentum between the 
protons, so that the correlated (C) and uncorrelated (Id) 
measurements of protons 1 and 2 are defined by 

C(Ap) = J P(p 1 ;p 1 + 2Ap) d Pl dfl p , 

U(Ap) = 1 J P ( Pl ) P ( Pl + 2Ap) d Pl dn p . (14) 

The integration in Q p is over all orientations of Ap. 
P (p) = / P (Pi P') dp' is the probability to measure the 
momentum p for one of the protons, irrespective of what 
the momentum of the other proton is. N is the total 
number of particles measured, i.e., N = j P (p) dp. 
The new correlation function is defined as 



C(A P ) 



C(A P ) 
U{Ap) 



1. 



(15) 



In figure [2] we show the correlation function for 45 Fe 
decay with E — 1.1 MeV, as a function of the relative mo- 
mentum of the protons. The short-dashed curve is for the 
singlet state and all others for the triplet state. The dot- 
ted curve does not include the proton-proton interaction, 
the long-dashed curve includes the Coulomb interaction 
between them and the solid curve includes both nuclear 
and Coulomb interaction between the protons. 

One sees that the properties of the correlation func- 
tions in the singlet and triplet states are completely dif- 
ferent. When Ap is small the correlation function is neg- 
ative only for the triplet state. It is -1 at Ap = for 
the triplet state, whereas it is close to zero for the singlet 
state. While for the former case the correlation function 
crosses zero at two points, it does not have a null point 
for the singlet case. It is also worthwhile mentioning 
that the effect of the Coulomb interaction between the 
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FIG. 2: Correlation function, eq. 1 151 for the 45 Fe decay with 
E — 1.1 MeV, as a function ot the relative momentum of the 
protons.. The short-dashed curve is for the singlet state and 
all others are for the triplet state. The dotted curve does not 
include the proton-proton interaction, the long-dashed curve 
includes the Coulomb interaction between the protons and 
the solid curve includes both nuclear and Coulomb interaction 
between the protons. 



protons can be switched off and the resulting correlation 
function C(Ap) multiplied by the Gamow factor in cq 
[13] yields a result (not shown in fig. [2]) slightly different 
than the long-dashed curve. In fact, the Gamow factor 
of eq. [T31 tends to underestimate the Coulomb final state 
interaction between the protons. 

The method described above is directly applicable to 
determine the spin mixing of final states in low-energy 
two-proton nuclear decay for + — > + transitions. In 
this case the final spin wave function of the pair equals 
that of the initial wave function. In particular, when 
singlet states are identified, spin-spin coincidence exper- 
iments will generate dichotomic outcomes for each single 
measurement. 

Fig. 2] shows the correlation function, C (Ap), for ro = 
4 fm and for different admixtures of singlet and triplet 
states. The dotted, dashed and solid lines correspond to 
M. = 0.1, 0.5 and 0.9, respectively. M is the absolute 
contribution of the triplet state. One clearly sees that 
different admixtures lead to very different dependence 
on Ap. 

Summarizing the above, we can say that the strong 
and Coulomb final state interactions cannot be neglected 
when the volume of spatial separation of the two-proton 
wavefunction is measured by ro < 6 fm. When ro < 4 
fm the strong final state interaction is noticeable in the 
relative motion spectrum of the two-protons (see fig. [I]) 
and its presence is reflected in the strong reduction of 
the relative energy. The tidal Coulomb force due to the 
charge of the daughter nucleus tends to increase consid- 
erably the relative motion of the two protons. The same 
applies for the Coulomb repulsion between the protons 
due to their own charge. These results point to the im- 
portance of considering FSI in the experimental analy- 
sis of two-proton decay experiments. The effect of final 
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FIG. 3: Correlation function, C (Ap), for two-proton triplet 
state decay of 45 Fe with ro = 4 fm (solid curve) and ro = 10 
fm (solid curve). 




10 20 30 

Ap [MeV/c] 



50 



FIG. 4: Correlation function, C (Ap), for ro = 4 fm and for 
different admixtures of singlet and triplet states. The dotted, 
dashed and solid lines correspond to M = 0.1, 0.5 and 0.9, 
respectively. M is the absolute contribution of the triplet 
state. 



state interactions are also visible in correlation functions 
which are considerably modified as the initial separation 
of the two-protons are probed. This is shown in fig. [31 
where the correlation function, C (Ap), is plotted for two- 
proton triplet state decay of 45 Fe with ro = 4 fm (solid 
curve) and ro = 10 fm (solid curve). 

We have shown that correlation functions, when de- 
fined appropriately, can clearly resolve the statistic na- 
ture of the diproton spin state in nuclear decay. This 
paves another route to study important problems of ba- 
sic quantum mechanics interest, such as the Einstein- 
Podolski-Rosen paradox [14j and Scrodinger cat states. 
Two-proton radioactivity may supply yet another test of 
the EPR dilema, namely whether the spin state function 
provides a complete description of quantum mechanics as 
we know it. These studies would be complementary to 
others performed in quantum optics and atomic physics, 
as well as in some instances of nuclear physics [161 Il8 |. 
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